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Explicit formulas for the solutions of the recurrence relations for 3-loop vacuum integrals are 
suggested. This formulas can be used for direct calculations and demonstrate a high efficiency. 
They also produce a new type of recurrence relations over the space-time dimension. 

- - - 1 Introduction 

Recurrence relations are powerful tools for evaluating multi-loop Feynman integrals [jj. They relate 
Feynman integrals with various degrees of their denominators. In many cases they provide the possibility 
to express an integral with given degrees of denominators as a linear combination of a few master integrals 
5— i | with some coefficients which we will call weight factors. 

At two-loop level the recurrence relations are relatively simple and one can easily find and realize the 
corresponding recursive algorithm for the calculation of the weight factors. 

At three-loop level the recurrence relations are more complicated and to find an algorithm to calculate 
weight factors is a serious problem. 

For vacuum integrals with one nonzero mass and various numbers of massless lines the corresponding 
algorithm was described in Pj-pJ- Here a repeated application of various recurrence relations is performed 
C*") ' until the integrals ofa certain type are eliminated. 

| In recent works [@],||-[7| the recursive algorithm was successfully applied in calculations of two-loop 

t^"- , and three-loop QED and QCD. 

\Q ■ Nevertheless, such recursive algorithms lead to too time and memory consuming calculations because 

of the size of intermediate expressions grows exponentially with respect to the degrees of the denominators 
in the initial integral. In fact, the calculations mentioned above were made at the limits of computer 

C capabilities. 

\Q ■ In this work we suggest a new approach based on explicit formulas for the solutions of the recurrence 

' relations. As an application, the case of three loop vacuum integrals with four equal mass and two massless 

lines is considered. The efficiency of this approach is demonstrated by calculations of previousely unknown 
coefficients in Taylor expansion of QED photon vacuum polarization for small q 2 . 



2 General case 



Let us consider the three-loop vacuum integrals with six different masses: 

m 2E i"*- 3I) f f f d D pd D kd D l 
B{n, D) = B(m, i%2, Ti3, Ti4, n>5, n&, D) = — -3 / / / nni n n 2 n „ 3 n „ 4 nns n » 6 (1) 



> 
• 1— 1 

X 

[lTT D / 2 r(3- D/2)] J J J D l D 2 D 3 D 4 D 5 D 6 

where 

D 1 = k 2 - ^m 2 , D 2 =l 2 -^2m 2 , D 3 = (p + k) 2 - fi 3 m 2 

D 4 = (p + I) 2 - [j, 4 m 2 , D 5 = (p + k + I) 2 - fi 5 m 2 , D 6 = p 2 - fi 6 m 2 
Let us derive recurrence relations that result from integration by parts, by letting (d/dpi) ■ Pj act on 
the integrand, with pi t j £ {p, k, I}. For example, acting by (d/dk) ■ (p + k) we get 



(D-2n 3 -n 1 -n 5 )}B(n,D) = {rnl + (3 -6 + ^ 3 - (i 6 + t n) + 2n 3 3+/i 3 

+n 5 5+(3- - 2 + M3 - + Ws)}-B(a, D ) ( 2 ) 

where l ± i3(ni, . . .) = B{ni ± 1, . . .), etc. 

Others relations can be obtained from (|^) by proper permutations of the n-j, /it and objects. 
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The common way of using these relations is step by step re-expression of the integral (gj) with some 
values of rij through a set of integrals with shifted values of n, , with the final goal to reduce this set to a 
few integrals with are equal to or 1, so called "master" integrals. The result can be represented as 

B(n,D) = J2f k (^D)N k (D) 

k 

where the index k enumerate master integrals Nk(D) and corresponding coefficient functions f k (rii,D). 

There are two problems on this way. First, there is no general approach to construction of such 
recursive procedure, that is to find proper combinations of these relations and a proper sequence of its 
use is the matter of art even for the case of one mass Q. Second, even in cases when such procedures were 
constructed, they lead to very time and memory consuming calculation because of large reproduction 
rate at every recursion step. For example, the relation (Q) expresses the integral through 7 others. 

Instead, let us construct the coefficient functions f k (n, D) directly as solutions of the given recurrence 
relations. 

For that, let us diagonalize the recurrence relations with respect to operators. We found that 

the recurrence relations can be represented in the following simple form 



{P(x 1: ...,x 6 )- n i l + — di(P(xi, . . ,,x 6 ))} x . =I - +fi .B(n,D) = 0, i = 1,...,6. (3) 

where 

P(xi,...,x 6 ) = 2(x 1 x 2 (x 1 + x 2 ) + x 3 x 4 {x 3 + x 4 ) + x 5 x 6 (x 5 + x 6 )) 
+xix 3 xe + X1X4X5 + X2X 3 x 5 + x 2 X4Xe 
— {x\ + x 2 + x 3 + X4 + £5 + xe)(xiX2 + x 3 X4 + X5Xe) 

The differential equation corresponding to (||) has the solution P D l 2 ~ 2 {x,i + fii). Let us consider 
"Lourent" coefficients of this function: 



1 I I I I I I dxidx 2 dx 3 dx4dx 5 dx 6 -.012-2 lA \ 

I{n - D) = Wftffftf xTxTxTx^xTx^ P{X1 + J*. •••■*« + /*) ( 4 ) 

where integral symbols denote six subsequent complex integrations with contours which will be described 
below. If one acts by (Q) on (^) one gets up to the surface terms the same expression as acting by 
Pdi -(D/2- 2){diP) on p D l^^ that is zero. Then, the surface terms can be removed if we choose 
closed or ended in infinity point contours. For more accuracy one can consider analytical continuations 
of the result on D from large negative values. So (Q) is the solution of the relations (g), and the different 
choices of the contours correspond to different independent solutions. Note, that if one chooses the 
contour as a small circle over zero, one get the true Lourent coefficient of the function F D I 2 ~ 2 , so this 
function can be called generalized generating function for the solutions of the relations @j. 

Then, in accordance with the dimensional regularization rules, the integrals (Oj) are not equal to zero 
only if at least three among rii are positive. So it is natural to construct the solutions from those that are 
equal to zero if the index from definite three-index set ("Taylor" indexes) is not positive. One can obtain 
such solutions if one chooses the contours, corresponding to these indexes, as circles over zero. In this 
case these three integrations can be evaluated and lead to coefficient in the common Taylor expansion in 
corresponding variables. 

The three remaining integrations in general case lead to the sum of generalized hypergeometric seria, 
but for some cases of practical interest (see below) can be reduced to the finite sums of Pochhammers 
symbols products. 

3 Example 

As an example let us consider the case of integrals with four equal mass and two massless lines, that is 
Mi = A'2 = 0, /is = /Z4 = (i5 = He = 1. Let us calculate the coefficient functions which corresponds to the 
choice of the master integrals from M. That is, we expand B(n) as 
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B(n, D) = N(n, D)B(0, 0, 1, 1, 1, 1, D) + M(n, D)B(1, 1, 0, 0, 1, 1, D) + T{n 7 D)B(0, 0, 0, 1, 1, 1, D) 

with the following normalization conditions 

iV(0,0, 1,1,1,1, -D) = 1, iV(l,l,0,0,l,l,-D) = 0, AT(0,0,0,1,1,1,£)) = 0, (5) 
M(0,0,1,1,1,1,-D) = 0, M(l, 1,0,0,1, 1,1?) = 1, M(0,0,0,1,1,1,D) = 0, (6) 
T(0,0,1,1,1,1,-D) = 0, T(l, 1,0,0,1,1,.D) = ; T(0, 0,0, 1, 1, 1,D) = 1, (7) 

The practical rule for choosing the integration contours is: circle around zero for unity in the master 
integral and contour over cut for zero in the master integral. 

To get N(n) one should make first the Taylor expansion in 2:3, £4, X5, xq 



/; '\/ ./ ^( (n3-l)!...K-l)! ^^^ 3 + 1. • • • + V D/2 ~ 2 ) 
The remaining integrals over xi, x 2 are of the type 



where we follow the normalization (ph. 

The case M(n, D) is analogous. The only difference is that due to the symmetry of the task we should 
take the sum of the solutions with the signatures (+ + ± ± ++) and (+ + + + ±±). 

The case T(n, D) is more complicated. The symmetry of the task assumes that one should try the 
following form of the solution 



T(ni, n 2 , n 3 , 774, n 5 , n 6 , D) = t{n\, n 2 , n 3 , 774, n 5 ,n 6 , D) + t(ni,n 2 , n 4 , n 3 , n 6 , n 5 , D) 

+t(n 1 ,n 2 ,n 5 ,n 6 , 773, 72,4, D) + t(n u n 2 , n 6 , n 5 , 774, n 3 , D) 

where t(n,D) is non-zero only if 774,775,776 > 0. Let us construct t(n,D) using (|]), keeping in mind 
possible mixing with N(n, D) solution. After differentiating over last three indexes the task reduces to 
the construction of t(n\, n 2 , TI3, 1, 1, 1, D). Let us consider the corresponding integral: 



-/ 1 / / / dX\dx 2 dx 3 , n . ,,, n/n_o , . 

t(ni, n 2 , n 3 , D) = ^— ^ <J> <J> <J> ^T^^i ( x 3 ~ xix 2 x 3 + x 1 x 2 (x 1 + x 2 - 4)) 1 (8) 

For 77,3 < 1 one can calculates this integral immediately (the possible N(n,D) contribution vanish). 
Taking into account the normalization (Q) we get 



f(n 1 ,n 2 ,n 3 <l,l,l,l,-D) = { - t(0 ^ 

(2-_D) ( „ 1+ „ 3) (2-_D) (t12+ „ 3) (-^-^-) ( _„ 1 _„ 3) (-^j^-) ( _„ 2 _„ 3) 
(-4)<™i+™2)(-8)™3 

[-ns/2] n-1 . ,D-1 . , . 

^ (3=^)_ fc (i)_»(-n s -2*)l 

For 773 > 1 using integration by parts for x 3 in (^) (which reduces to evaluation of (773 — \ ) th derivative 
Of P^/ 2 - 2 ) the t(77l,77 2 ,773, D) can be reduced to a set of t(ni, 772, 1, D) with different ni,n 2 . Let us 
extract the t(ni, n 2 , 1, 1, 1, 1, D) from ?(m, 772, 1, D) according to the conditions (0) 



t( ni ,n 2 , 1, 1, 1, 1, D) = j ^ (i(n 1 ,n a , 1, D) - 1(0, 0, 1, D^m.na, 1, 1, 1, 1)) (9) 
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One can calculate the t(ni, n 2 > 1, 1, 1, 1, D) by direct use of the (g) expanding it for example in seria 
over D/2 — 2, but we found more suitable to use the recurrence relations on ni,n 2 : 

t(ni,n 2l l, 1,1,1,£>) = - 4 (J3-lj( 2 n 2 1 ) --D+2) ~ 1^2 - 1, 0, 1, 1, 1, £> - 2) 

+t(m- 2, 7*2-1, 1,1, 1,1,^-2)) 

2(L>-2) 2 (ll_D-38) AT ( i i i i /, n \ 

~ 3(3g-10)(3£>-8)(£>-3) N ("1 , "2, 1, 1, 1, 1, -D) (10) 

t(n x , n 2 ,l,l,l,l,£>) = Sr-^j tK, n 2 + 1, 0, 1, 1, 1, D) 

+ {ifE^ (ni-l,n 2 + l,l,l,l,l,D) (11) 

With the help of (EJ) the ni + rt2 can be reduced to — 1, 0, 1 and with the help of ([ll]) the m — ri2 can 
be reduced to 0, 1 (noTe that t(ni, n 2 , 1, 1, 1, 1, D) = t(i%2,ni, 1, 1, 1, 1, D)). Here at every recursion step 
the one integral reexpreses through the other one plus rational over D, that is there is no "exponential 
reproduction". Then, the recursion acts separately on variables n\ + n-i and n\ — n 2 . So, although the 
relations (0,0) can be solved to explicit formulas, this "safe" variant of recursion is in this case the 
most effective way of calculations. 

The relations ( JTo| , |TT| ) are the simple example of the recurrence relations with D-shifts, which can be 
derived in the following way. Note that if f k (n t ,D) is a solution of (g), then P(I~ + Hi)f k (n t ,D - 2) 
also is a solution. Hence, if f k (rii, D) is a complete set of solutions, then 



f k (m, D) = S k n (D)P(\- + m)n?H,D - 2) (12) 



where the coefficients of mixing matrix S depends only over D. For the solutions (Q) the matrix S is unit 
matrix. On the other hand, the desire to come to some specific set of master integrals leads to nontrivial 
mixing matrix and for the example considered above these coefficients are 

an _ 3 (3D-8)(3g-10) gm _ 3 (3D-8)(3D-10) qt _ (D-2) 2 

°n 64 (D-4) 2 °rn 16 (2D-7)(2D-9) °t 4(D-3)(£>-4) 

qt (11-D — 38)(,D — 2) ™ qt qm qm qn r\ 

J n — 32(D-3)(D-4) 2 a t ~ °m — °t ~~ J n ~ °m ~ u 

To check the efficiency of this approach we evaluated, to 3 loops, the first 5 moments in the z = 
q 2 /4m 2 — > expansion of the QED photon vacuum polarization 



n(z) = ^C„z" + 0(a 4 ), 



n>0 



The C n are expressed through approximately 10 5 scalar integrals, but there is no necessary to evaluate 
these integrals separately. Instead, we evaluated a few integrals of (|j) type, but with p D /' 2 ~ 2 producted 
by a long polinomial in xi. 

After OS mass || ^| and chargeQ renormalization, we obtained the finite D — > 4 limits (the coefficients 
Ci , C 2 , C3 can be found in Q ) : 

C //V 2 r256/- I 2522821 f 129586264289 \] 

^ 4 \ ly L 693 S2 ^ 9437184 S3 143327232000 

3 

, at [160 (-1 8 i q\ f , 1507351507033 f 269240669884818833 1 1 a 
~^ L 231 V 5 / ^ 2 1651507200 ^ 3 245806202880000 J / ^3 

2 

_j_ 51986 N ^_ , _32_ ju« 



127575 ^2 ' 693 _ 



/-( _ f at 2 [ 1024 a , 1239683 a 512847330943 \] 

^ 5 \ y L 3003 ^ 2 3932160 ^ 3 556351488000 ''J 

, at [ 640 /-i 8 l o\ /■ I 939939943788973 f 360248170450504167133 1 1 a 

" r LlOOl V 5 / 190749081600 S3 60837035212800000 J J 7] _. 



a 2 



, 432385216 AT_ i 512 iy~ 

^ 1260653625 n 2 ~*~ 15015 _ 
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where we follow common practice 10 , by allowing for N degenerate leptons. In pure QED, N = 1; 
formally, the powers of N serve to count the number of electron loops. 

The N contribution of C4 is in agreement with recent QCD calculations jjj, the N 2 part of C4 and 
the C5 are new. 

The bare (non-renormalized) integrals were calculated for arbitrary D. Calculations for C4 were made 
on PC with Pentium-75 processor by REDUCE with 24Mbyte memory, within approximately 10 CPU 
hours. The most difficult diagrams for C5 were calculated on HP735 workstation. 

These results demonstrates a reasonable progress in comparison with common recursive approach. 
For example, the common way used in Q demands several CPU hours on DEC-Alpha workstation to 
calculate full D dependence of Ci integrals, and further calculations became possible only after truncation 
in (D/2 — 2). In the present approach the full D calculations for C2 demand about 5 minutes on PC. 



4 Conclusions 

The new approach suggested in this work allows to produce explicit formulas (H) for the solutions of 
the recurrence relations for 3-loop vacuum integrals. This formulas can be used for direct calculations 
and demonstrate a high efficiency. On the other hand, they produce a new type £>-shifted recurrence 
relations (|l2|) for these integrals. Finally, we hope that simple representation (||) of the traditional 
recurrence relations which allows to obtain all these results is not intrinsic for 3-loop vacuum case and 
generalization for multi-loop or/and non- vacuum case is possible. 
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